o 
o 

(N 



< 
(N 



X 



SLAC-PUB-11400 
UMN-D-05-2 
SMUHEP/03-14 

Two-boson truncation of 
Pauli— Villars-regulated Yukawa theory^ 



Stanley J. Brodsky, John R. Hiller,^ and Gary McCartor*^ 



W). '^Stanford Linear Accelerator Center, Stanford University, Stanford, California 

^ ■ 94309 



^Department of Physics, University of Minnesota-Duluth, Duluth, Minnesota 

55812 

'^Department of Physics, Southern Methodist University, Dallas, TX 15215 



> 

in 

0^ 
(N 
00 . 

O . Abstract 

vn 

. We apply light-front quantization, Pauli- Villars regularization, and numerical tech- 

1 ^ ' niques to the nonperturbative solution of the dressed-fermion problem in Yukawa 

I ' theory in 3 + 1 dimensions. The solution is developed as a Fock-state expansion 

truncated to include at most one fermion and two bosons. The basis includes a 
r~| , negative-metric heavy boson and a negative-metric heavy fermion in order to pro- 

vide the necessary cancellations of ultraviolet divergences. The integral equations for 
the Fock-state wave functions are solved by reducing them to effective one-boson- 
one-fermion equations for eigenstates with = 1/2. The equations are converted 
I to a matrix equation with a specially tuned quadrature scheme, and the lowest 

mass state is obtained by diagonalization. Various properties of the dressed-fermion 
state are then computed from the nonperturbative light-front wave functions. This 
work is a major step in our development of Pauli- Villars regularization for the non- 
perturbative solution of four-dimensional field theories and represents a significant 
advance in the numerical accuracy of such solutions. 
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1 Introduction 



One of the important unsolved problems in obtaining bound-state solutions 
in quantum field theories such as quantum chromodynamcs is how to im- 
plement ultraviolet renormalization nonperturbatively. The central difficulty 
is that any truncation or approximation of the theory which breaks Lorentz 
symmetries, such as the truncation of the Fock space, introduces spurious di- 
vergences. Quadratic divergences can occur in the approximated theory, even 
when the renormalizable perturbative Feynman theory has only logarithmic 
divergences. Such problems arise even when one uses light-front quantization 
methods which retain the maximal set of kinematical Lorentz symmetries. 
A systematic light-front approach to renormalization, developed by Glazek, 
Harindranath, Perry, and Wilson [1,2,3], is to introduce new effective inter- 
actions which act as counterterms to control the ultraviolet behavior of the 
approximated theory. However, this method is challenging to implement in 
practice, because the large set of new effective interactions, some of which 
have a nonlocal structure, can require so much input data that the ability to 
make predictions is seriously compromised. 

Another important nonperturbative approach is to cast the bound-state prob- 
lem in the form of an effective Bethe-Salpeter equation and use Schwinger- 
Dyson methods to construct the effective renormalized quark and gluon prop- 
agators of the theory, including running quark masses consistent with chiral 
symmetry breaking [4,5]. In principle, this method could be used to predict 
the light-front wave functions needed for QCD phenomenology, if one evalu- 
ates the Minkowski-space Bethe-Salpeter wave functions at fixed light-front 
time. However, it has been difficult to carry out this program in practice since 
the analyses have only been done in ladder gluon-exchange approximation and 
in Euclidean space. The full structure of nonperturbative renormalization will 
require consideration of higher-order kernels. 

Lattice gauge theory [6] is currently the most successful method for solving 
gauge theories; it has a systematic gauge-invariant regularization procedure 
despite the fact that nonlinear gauge interactions are introduced and the lat- 
tice structure itself violates Lorentz symmetries at finite lattice size. How- 
ever, unlike light-front methods, wave function information for bound states 
and other physical features of hadron dynamics have to be obtained indirectly 
through moments, because the lattice theory is formulated in Euclidean space. 

In our work [7,8,9,10,11,12] on light-front quantization, we have shown that 
the technique of Pauli-Villars (PV) regularization [13] of ultraviolet diver- 
gences can be implemented in nonperturbative light-front calculations of field- 
theoretic bound states. Until now this has been limited to light-front Fock- 
space truncations where the resulting wave functions could be computed ana- 
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lytically [11,12] or to numerical calculations using discretized light-cone quan- 
tization (DLCQ) [14,15] with limited precision [9]. In this paper we will test 
the PV method in some detail and present accurate results for the dressed- 
fermion state in Yukawa theory with a Fock basis which includes three-particle 
states. 

In order to perform our calculations using PV regularization, we first intro- 
duce a sufficient number of PV fields in the underlying Lagrangian in order 
to render perturbation theory finite. However, we must also make sure that 
our nonperturbative result, if expanded in a power series in the coupling con- 
stant, would give agreement with the usual Feynman series for processes that 
can be calculated perturbativcly. Obtaining such assurance may require the 
introduction of additional PV fields or counterterms, or both. Paston and 
Franke [16] have given a general set of rules which determine exactly what 
combination of PV fields and counterterms are needed to assure perturbative 
equivalence with Feynman methods. In that paper they apply their methods 
to the case of Yukawa theory; in [17], Paston, Franke and Prokhvatilov apply 
the methods to give the same information for QCD. In the case of the Yukawa 
theory that we study here, one PV fermion and one PV boson are sufficient 
to assure perturbative equivalence with Feynman methods without the need 
for counterterms. The PV regulators always preserve Lorentz invariance and 
may preserve gauge invariance, or at least the Ward identity. In the theories 
studied to date, the PV fields alone have been sufficient to give a finite the- 
ory, perturbativcly equivalent to Feynman theory, without breaking Lorentz 
or gauge invariance; in future work we expect to include counterterms, along 
with the PV fields, in order to restore symmetries broken by the PV fields. 

We must truncate the Fock space in order to be able to perform a numerical 
calculation of the spectrum and the wave functions. So far, all of our calcula- 
tions have been for the ground state. The truncation breaks the symmetries 
of the theory, but since this introduces only finite breaking, we argue that if 
our answer is close, numerically, to the answer without truncation, then the 
approximated result should yield a useful result. This was shown to be the 
case in a nonperturbative calculation of the electron's magnetic moment in 
QED [12]. 

We will work with the theory in light-cone coordinates [18], in order to have 
a well defined Fock expansion and a simple vacuum [15]. The field-theoretic 
bound-state problem for the dressed fermion is then reduced to integral equa- 
tions for the wave functions which appear in the Fock expansion. The ex- 
pansion is truncated to include states in three sectors: the bare fermion, the 
one-boson-one-fermion states, and the two-boson-one-fermion states. We ne- 
glect fermion-pair contributions. 

In Yukawa theory and in QED, both without fermion loops, it is adequate to 
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include one PV fermion and one PV boson to regulate the ultraviolet diver- 
gences. At the end of the calculation, we wish to take the PV masses large, 
and there is an ambiguity as to how to do that. We can take the fermion PV 
mass much larger than the PV boson mass, take the PV masses to be about 
equal, or take the PV boson mass to be much larger than the PV fermion 
mass. The answer that one obtains depends on the choice of this mass ra- 
tio [11]. In the case of QED, restoration of gauge invariance requires that the 
PV fermion mass be taken to infinity while the PV boson mass remains fi- 
nite [12]. (Although both types of PV particles must be included to render the 
theory finite, once the calculations are complete the limit of the PV fermion 
mass going to infinity exists.) Of course, since the limit of the PV fermion 
mass going to infinity is finite, taking the fermion mass sufficiently large but 
finite will be adequate. We believe that taking the fermion mass to infinity 
first is the physical limit in Yukawa theory as well, but we have no argument 
as strong as the one in the case of QED; in the present paper we shall show 
results for taking the PV fermion mass large first and also for keeping the PV 
masses equal. Due to the truncation of the representation space, there always 
remain uncancelled divergences, so at least one PV mass must remain finite. 

Unlike QED, where the wave function and vertex renormalizations cancel, 
Yukawa theory has a logarithmically divergent charge renormalization even 
without considering fermion loops and vacuum polarization. Thus in our trun- 
cation we must consider renormalization of both charge and fermion mass. To 
handle these renormalizations nonperturbatively, we regulate the theory and 
then fix the bare coupling and bare fermion mass by imposing conditions on 
the mass and Dirac radius of the dressed-fermion eigenstate. 

There are three objectives for the present work: One is associated with the 
need to take the limit of the PV fermion mass to infinity first, leaving the 
theory regulated by just the PV boson mass. While the answer for infinite PV 
fermion mass exists, an explicit solution is not known, and for now we must 
settle for taking the fermion mass very large compared with the PV boson 
mass. That choice of masses raises a strong possibility of numerical trouble. 
Numerical methods work best in problems with a single scale, but here we 
have three scales: the physical mass scale, the PV boson mass scale, and the 
PV fermion mass scale. In this paper we present techniques which allow us to 
make accurate calculations even when the three mass scales are very different. 

A further objective of the present work is associated with the need to keep 
one PV mass, usually the PV boson mass, finite, and therefore the need to 
choose a value for it. In [12] we proposed a method for choosing the PV boson 
mass based on the idea that there are two types of error associated with a 
finite PV mass, and that we should choose a value for which neither type of 
error is too large. (If no such value exists then our method will not work.) Two 
elements of the proposed method of choosing a final PV mass are examined in 
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the present paper. One proposal is that one can obtain a reasonable estimate 
of the weight of the true wave function on the lowest excluded Fock sector 
by performing a perturbation calculation using the calculated wave function 
as the unperturbed state. The other proposal is that the weight of the wave 
function on the first excluded sector provides a useful estimate of the percent 
by which physical parameters would shift if the first excluded sector were to be 
included. We use the results of a calculation including only up to two particles 
to perform a perturbative estimate of the weight of the true wave function 
on the three-particle states and then compare that result with the weight of 
the wave function in the three-particle sector when that sector is included in 
the calculation. Wc also calculate a number of physical quantities for both 
the two-particle truncation and the three-particle truncation and compare the 
percent by which they change with the percent of the wave function that is in 
the three-particle sector. 

The final objective of the present work is just to study the effects of adding 
an additional sector to the calculation. This extends earlier work [11] where 
a truncation to two particles allowed analytic solutions. We can see whether 
properties of the solutions persist with inclusion of the three-particle sectors 
and whether new physics appears. 

The PV fermion and PV boson are introduced to the theory in such a way 
that the interaction term in the modified Lagrangian is a coupling between 
zero-norm combinations of the physical and PV fields. This guarantees that 
instantaneous-fermion terms do not appear in the light-cone Hamiltonian and 
only ordinary three-particle vertices remain. The negative-metric particles also 
cancel infinities from integrals over transverse momenta. The numerical ap- 
proximation is then applied to a finite theory at fixed PV masses, and the be- 
havior of the solution is studied as the PV masses are increased. The structure 
of the problem is simpler than in the case of quantum electrodynamics [12], 
not only because the bosons are massive scalars but also because the kernels 
of the integral operators are not plagued by singularities caused by spurious 
thresholds. Here the bare fermion mass is driven to large values rather than 
the small values seen in QED; this prevents the singularities. 

The solution to the eigenvalue problem for the dressed-fermion state is ob- 
tained by first rearranging the coupled set of integral equations for the mass 
eigenvalue problem into a set of effective equations for the one-boson-one- 
fermion wave functions. The kernels of these equations have contributions from 
intermediate states containing a single bare fermion and states containing two 
bosons and a fermion, including self-energy terms. The effective equations can 
be considered an eigenvalue problem for the bare coupling at fixed bare and 
dressed fermion masses. The numerical solution provides the wave functions 
as well as the bare coupling. We can use these to compute various properties 
of the dressed fermion, including the magnetic moment and axial coupling, as 
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well as structure functions. 



The traditional DLCQ method [14,15] has limited accuracy in the present 
case, because the solution is sensitive to regions of small longitudinal light- 
cone momentum fractions, on the order of the reciprocal of the PV mass 
squared in units of the physical mass. The equal spacings in momentum used 
in DLCQ must then be so large in number as to be impractical. Here we use 
coordinate transformations and Gauss-Legendre quadrature to capture these 
important regions. 

The calculation of the solution is simplified by working in transverse polar 
coordinates. The Fock expansion is constructed to be an explicit eigenstate of 
Jz- The dependence of the wave functions on the azimuthal angle can then be 
determined exactly and removed from the numerical calculation. This reduces 
the effective dimension of the numerical problem from three to two. 

The notation that we use for light-cone coordinates is 



The time coordinate is x"^, and the dot product of two four- vectors is 



The light-cone momentum component conjugate to x~ is p"*", and the light- 
cone energy is p~ . Light-cone three- vectors are identified by underscores, such 
as 



For additional details, see Appendix A of Ref. [7] or the review [15]. 

We begin in Sec. 2 by introducing the light-cone Hamiltonian and dressed- 
fermion Fock-state expansion for Yukawa theory and by giving expressions for 
quantities to be computed from the Fock-state wave functions. Previous results 
for the one-boson truncation [11] are summarized in Sec. 3. Section 4 contains 
the analysis of the two-boson truncation, including a summary of the numerical 
techniques and a presentation of the results. Some concluding remarks are 
given in Sec. 5. Details of the kernels and the numerical approximation are 
left to two appendices. 

Somewhat related calculations have been done by Bylev, Glazek, and Przcs- 
zowski [19], except that they did not use a covariant regulation procedure. 
Similar work in a purely scalar theory has been done by Bernard et al. [20], 



x"^ — x'^ + x^, x± — {x^,x'^). 



(1) 




(2) 



(3) 



6 



and more recently the two-fermion problem has been considered by Mangin- 
Brinet et al. [21]. For a more formal treatment of dressed constituents, see [22]. 



Other light-cone methods that show promise include supersymmetric DLCQ 
(SDLCQ) [23] and the transverse lattice approximation [24]. Both are used 
specifically for gauge theories. 



2 Yukawa theory 



We consider Yukawa theory with a PV scalar and a PV fermion. The action 
is 



(4) 



lYd,, - (d^i^ih^) ^1 + "^i^iV'i - ^(00 + 0i)(^o + i'i)i^o + ^i) 



The subscript indicates physical fields and 1, PV fields. The fermion masses 
are denoted by rrii, and the boson masses by nj. When antifermions are ex- 
cluded, the resulting light-cone Hamiltonian is 



i,s •' P 



P 

,'2 I 



(5) 



+ H / dP^l { [^-2siP^ g) + V2s{p + q, q)] b]Ap)al{q)bi_s{p + g) 
+ [Uj(p,q) + Ui(p + q,q)] + £) + /i.e.} , 



where a\ creates a boson of type j, bj ^ creates a fermion of type i and spin s, 



Uj{p,q) 



9 ^3 9 ^2s ■ P± 



Vl67r3p+v^ 



-:, V2s{P,q) 



(6) 



and 



e2s 



'V2 



(2s, i). 



(7) 
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The nonzero (anti)commutators are 



h(g),a](£')]=(-irM(^-g'), (8) 

The opposite signature of the PV fields is the reason that no instantaneous- 
fermion terms appear in P~ ; these terms are individually independent of the 
fermion mass and cancel between instantaneous physical and PV fermions. 

We expand the eigenfunction for the dressed-fermion state in a Fock basis as 



$+(£) = E ^iblim) + E / dqfijsiq)bUP - q)a]iq)\0) (9) 

i ijs 

dqidq^fijks (gi , ^2 ) , -= h\s (Z - gi - 52) a] (^i) 4(^2)10) + . . . 



and normahze it according to $^ ■ $0- = 6{P^ — P). The wave functions / that 
define this state must satisfy the coupled system of equations that results from 
the field-theoretic mass-squared eigenvalue problem P~^P~^^ = M^$+, since 
we work in the frame where P± is zero. The state has — +1/2. The first 
three coupled equations are 



p+ 

i\z, + E(-l)"+^i^+ / dq {U,,.{q){VAP -q,q) + V1{P, q)] 

+Uj+{q)[U,{P-q,q) + U,{P,q)]] = Mh„ (10) 



1 -y y 



ks{q) + E(-l)^' {^^'K-[y;{P -q,q) + V-{P,q)] 

i' 

+Zi'SsAUi{P -Q:Q) + Ui'{P,q)]} (11) 
+2EV=^^^ / c?g'{^,fe,_.(g,g')[^2.(P-£-£',£') 

i',k Y 1 + ^jk 



and 
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l-yi-y2 yi 1/2 



fijks{Ql,Q2) 



(12) 



i' 



X [K!2,(Z -01-02,22) + V2s{P - gi, 02)] 

+ fi' j siOiWiiE - qi- 02,02) + Ui>{P- 01,02)] 

+fi'k,-s{g2)[Vl2siP -01-02,01) + V2s{P - gi)] 
+Uks {02) [U^ {P - Ol- 02,01) + Ui'iE- 02,01)]} + ■■ ■ 

= M^fijksiqi,q2)- 



The equations are invariant under Lorenta boosts: P~^ 7P"'". We represent 
these diagrammatically in Fig. 1. 

Each wave function has a total eigenvalue of for s = +1/2 and 1 for s = 
— 1/2. For the one-boson wave functions, this greatly restricts the dependence 
on the azimuthal angle; however, for the two-boson wave functions, the total 
Lz eigenvalue can be obtained in an infinite number of ways by combining 
different individual orbit als. 



Fig. 1. Diagrammatic representation of the first three coupled equations for the 
wave functions of the dressed-fermion state, Eqs. (10)-(12) in the text. The large 
blobs represent wave functions, the crosses represent light-cone energies, the solid 
lines indicate the fermion constituent, and the dashed lines correspond to bosons. 

We define the physical wave functions as the coefficients of Fock states contain- 
ing only particles with positive-definite norm. This reduction can be achieved 
by requiring all Fock states to be expressed in terms of the positive-norm 
creation operators b^g and aj and the zero-norm combinations bl = + b\g 
and = Oq -|- a|. We then discard any term containing a bl or an a^, which 
would be annihilated by the PV-generalized electromagnetic current, leaving 
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the physical state 



+ i:Jdq (^i-iy^' hs{q?j bUP-q)al{q)\0) 

f ( (-iy+3+^ \ 

+ 12 dqidq2 / fijksiqi,q2) 

s •' \ijk a/1 + "ifc / 



(13) 



X bl{P -qi- g2)4(?i)«o(?2)|0) + . . . 



It is normalized by 



+ 12 j dqidq2 



E(-i) 



i+j+k_ 



V2 



--fijksiOl, Q2) 



+ ... 



ijk \J^ + ^jk 

Prom this state we can compute a boson structure function 



(14) 



fBs{y) = J dq6{y-q+/P+) 
+ f dq,dq^Y.5{y-ql/P^ 

•' n=l 



E(-i)^^v.,.(?) 



1-3 



(15) 



E(-ir 

ijk 



1 + 57fc 



fijksiOl, 02) 



+ 



which is defined as the probability density for finding a constituent boson of 
longitudinal momentum fraction y when the constituent fermion has helicity 
s. The form factor slope F[{0) is given by [7] 



s ■' 

~E / dQidqiYl 



|v^E(-i)^^V.,.(i) 



(16) 



yi 



VuE(-i) 



i+j+k 



x/2 



ijk 



1 + 6. 



--fijksiQl, ^2) 



with Vi± = Xt^ — \-yir~- Since no transverse cutoff is used, this expression is no 
longer the approximation that it was in earlier work [7,8,9]. The Dirac radius 
of the state is given by = y^— 6F{(0). Prom the two-body wave function, we 
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compute a distribution function 



(17) 



and its moments 



J. J. 
/m = y dx(f){x,Q±), x = j dxx(f){x,Q±)/ fM, 



We can also calculate the axial coupling constant 



(18) 



gA = {zo-z,f + J2{-ir-'/'Jdq 



E(-i) 



i+j+k _ 



V2 



ijk ^JY+~5^ 
and the anomalous magnetic moment [25] 



fijks{qi,q2) 



+ 



(19) 



K = -M 



-M 



I 



yi 







\ijk 






^9qiyJ 


\dqix 



1 + 5jk 



(20) 



^ fljksiOi^Ql) 



^1 + ^jk 



Here an up or down arrow indicates wave functions associated with a dressed 
fermion having a value that is positive or negative, respectively. The wave 
functions with an up arrow are the same as those without an arrow; the coupled 
system of equations is constructed and solved for this case. The wave functions 
for the opposite spin are related as /i± = l^fl%, as can be estabhshed by 
comparing the equations that they each satisfy. 



Once the coupled equations are truncated to a finite system, we have a well- 
defined problem with cancellations of infinities between any infinite integrals. 
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The PV particles are kept in the basis to provide these cancellations. A one- 
boson truncation produces an analytically solvable problem, which we explored 
in [11] and discuss briefly here in Sec. 3. Less severe truncations produce larger 
coupled systems that require numerical techniques for their solution. Given an 
accurate discretization and the consequent finite matrix eigenvalue problem, 
one can compute mass eigenvalues and associated wave functions. The bare 
parameters, i.e. the bare coupling g and the bare mass mo of the positive- 
norm fermion, can be fixed by fitting "physical" constraints. Here we specify 
the dressed fermion's mass M and radius R; the fit to a chosen value of R is 
determined by an iterative root-finding scheme. 

The discretization must be chosen carefully. When the PV masses are large, 
the integrals become sensitive to small momentum fractions, of order /Xq/ /^^ 
and nl/ml. This makes traditional DLCQ [14] impractical, because it di- 
vides the longitudinal momentum into segments of equal size; the number of 
segments then grows as fJ^I/fJ^l and ml/fil, making the matrix diagonaliza- 
tion problem impossibly large. We instead employ a discretization based on 
Gauss-Legendre quadrature and certain variable transformations, as discussed 
more fully in Sec. 4 and Appendix B. The derivatives needed for computing 
the radius and the anomalous moment are estimated from cubic-spline fits 
to (m^ -|- q±)/{l — y) + (A*j + <fi_)/y fijs{y, Q±) as a function of q± at fixed 
y. The use of the multiplier reduces the variation in / and makes possible a 
better fit. Standard finite-difference approximations are not useful, because 
the quadrature scheme introduces very unequal spacings. 



3 One-boson truncation 



As an alternative to traditional DLCQ, we can explicitly truncate the system 
with respect to the total number of bosons in any Fock state. We have al- 
ready studied the case of the one-boson truncation [11] where the system of 
two equations can be solved analytically. The structure of the solution is as 
follows. From the second equation in the system, Eq. (11), the one-boson wave 
functions are immediately found to be 



P+ 



M2 



l-q+/P+ q+JP^ 



fij-iQ) 



P+ 

1-9+/P+ q+/P+ 



+ ^{-l)^z,UkiP,q) 

k 

I^Y.{-irz,'^V;{P-q,q) 



(21) 
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Substitution into the first equation, Eq. (10), yields algebraic equations for 
the bare-fermion amplitudes, which are 



(M^ - mf)zi^g'^iJ,l{zo - zi)J + g'^mi{zomo - zimi)Io 

+g'^IJ,o[{zo - zi)mi + zoruo - Zimi\li, (22) 



with 



167r' 1r M2 - ^^+'^^ - 1/(1 
1-2/ y 

i-y y 



1^0 



The presence of the PV regulators makes these integrals finite and allows Iq 
and J to satisfy the identity i^lJ{M^) = M^Io{M^). Because M is held fixed, 
these equations can be viewed as an eigenvalue problem for g^. The solution 
to this eigenvalue problem is 



2 (M =F mo)(M =F mi) zi M ^ mo 



^ ~~(mi-mo)(/io/i±M/o)' ~ M^mi' ^ 



Structure functions and distribution amplitudes can then be computed and 
the PV-mass limits studied. Figures 2, 3, and 4 show typical results. The 
two-boson contribution is computed perturbatively. The forms show a signifi- 
cant sensitivity to the Pauli-Villars masses. Further results can be found in 
Ref. [11]. 



4 Two-boson truncation 



4-1 Effective equations 



For a two-boson truncation, the solution is no longer analytic, but the coupled 
equations (10)-(12) can be reduced to eight equations for the two-particle 
amplitudes only, which are of the form 



1 -y y 



fijsiy, q±) = Yg-^ 2^ —^-——fajsiy, q±) 
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(a) 



(b) 
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Fig. 2. Bosonic structure functions fB±{y) as defined in Eq. (15) of the text, where 
y = xb = q'^/P'^, for the dressed-fermion state with mass M = jiq and radius 
R = 0.01///0- The wave functions are computed from a truncation to one constituent 
boson, with the two-boson contribution then computed pcrturbatively. The sohd hue 
includes both contributions. The long dashes show the one-boson contribution, and 
the short dashes show the two-boson contribution. For (a) and (b), the PV masses 
are mi = /xi = 2000;Uo, and for (c) and (d), they are mi = SOOOO^uq and /j^i = 500/xo. 

1 



167r2 



(26) 



dy'd(i'l4jlabAy^ v'^ (i'±)fabs'{y', q'±) 

-y 

Z / dy'dq'l4Jl^,,^,{y,qr,y',q'j_)fabs'{y',q'±), 
«' 

with the angular dependence removed via \/P^fij+{q) = fij+{y,q±) and 
VP~^fij-{q) = fij-{y, q±)e^'f'. Here lija is a computable self-energy and J,'-"'' 



9 1-y 



ijs,abs' 



is the kernel due to n-boson intermediate states. Specific forms for I and J^"^ 
can be found in Appendix A. A diagrammatic representation is given in Fig. 5. 



There is a difficulty with the structure of this eigenvalue problem, because the 
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(a) 



(b) 



Fig. 3. Same as Fig. 2 but for the distribution amplitude (j){x, Q±) defined in Eq. (17) 
of the text. The individual lines correspond to different values of Q±. For (a) the PV 
masses are mi = /xi = 2000/xo, and for (b) they are mi = 50000/xo and m = 500/xo. 
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Fig. 4. Moments of the distribution function, as defined in Eq. (18) of the text. 
The filled circles correspond to /m, the triangles to x, and the squares to x^. For 
(a) and (c), the PV masses are m\ = jii = 2000/xo, and for (b) and (d), they are 
mi = 50000/io and fj,i = 500/Lto. 
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Fig. 5. Same as Fig. 1 but for the effective equations in (26) of the text. 

original eigenvalue appears in the denominator of integrands that go into 
finding the matrix elements. If attacked directly, this nonlinearity requires 
an additional layer of numerical effort, to get a self-consistent solution. An 
indirect approach is to again convert the problem to one where is the 
eigenvalue. To maintain symmetry (up to the indefinite norm), this conversion 
is done by defining a new function 



+ ql ^ + gj _ ^2 



y 



(27) 



before completing the rearrangement. The smallest real g that is obtained for 
a given M is taken to be the coupling value for which M is the mass of the 
lowest state. Obviously, this works well with the renormalization condition 
where M is fixed as input. 



4-2 Method of solution 



We solve the reduced integral equations (26) numerically by converting them 
to a discrete matrix equation via Gauss-Legendre quadrature of the integrals. 
The resolution of the quadrature is characterized by the order of the underlying 
polynomial, which we denote by K in the longitudinal direction and by 2A^ + 1 
in the transverse direction. For the transverse quadrature, only odd orders are 
used, to keep gj_ = as a quadrature point, and only + 1 points are used. 
(The other N points correspond to negative qj_.) Thus K and N characterize 
the resolution of the approximation, and we consider large values in order to 
be close to the continuum limit. 

Before applying the quadrature rules, the variables y' and q'j_ are transformed 
in such a way as to emphasize those regions most important to the approxi- 
mation. In the case of q'j_, the transformation is chosen also to produce a finite 
range of integration without introduction of a cutoff. Since the PV contribu- 
tions make the integrals finite, no cutoff is needed for the continuum problem, 
and any cutoff would only be an artifact of the numerical approximation if 
introduced. A more complete discussion of the transformations is given in 
Appendix B. 
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The matrix eigenvalue problem is solved by applying the Lanczos diagonaliza- 
tion scheme developed previously [9]. This particular scheme was designed to 
efficiently handle the present situation where the matrix is self-adjoint with 
respect to an indefinite norm. What is different here, compared to the case 
in [9], is that the matrix is not sparse. Nevertheless, the Lanczos approach is 
much faster than standard diagonalization algorithms for nonsymmetric ma- 
trices, because we are interested in only one eigenstate. Typical matrix sizes 
are on the order of 20,000 by 20,000, but provide as many nonzero entries as 
the much larger sparse matrices considered in [9]. 

Use of the Lanczos technique is important for another reason. Standard di- 
agonalization routines have difficulty when there are multiple mass scales in 
the problem; specifically, when /Xi is intermediate between /iq and mi, such 
that fj,o <^ <^ mi, the standard diagonalization can fail. For the Lanczos 
process, we need only an accurate representation of the product of the matrix 
and a vector. The contribution of J^'^^ to the matrix can then be written in 
a factorizcd form Y^Ii^q^—^Y Vitvf,rj / {M"^ — where r/ is a diagonal ma- 

trix that represents the signature of the norm, as determined by the factors 
(—1)"+^ given as part of the definition of J^^"^ in Eq. (29). This factorized form 
provides a more faithful numerical representation of the contribution from J^^"^ 
by restricting the J^^^-vector product to a linear combination of the vectors 

Vi'. 



4-3 Results 

The convergence with resolution is indicated in Fig. 6. For most quantities, 
and for the range of PV masses considered, resolutions of = 50 and = 30 
are sufficient for convergence, and we use these resolutions for most of our 
calculations. 

The variation of the structure functions and some characteristic quantities 
with respect to PV mass can be seen in Figs. 7, 8, and 9. When the PV 
masses are equal, the radius of the state with mass M — /iq is driven to 
zero as the PV mass is made large; thus in this limiting case a fixed radius 
cannot be maintained. The distribution amplitude and its moments are shown 
in Figs. 10 and 11 for two sets of PV masses. 

A comparison of results for the one and two-boson truncations at fixed reso- 
lution and fixed PV masses is given in Figs. 12 and 13. The dressed-fermion 
radius R is held fixed at the same value of COl/^Uo for both truncations. 
The one and two-boson contributions to the structure functions are shown in 
Fig. 14; these should be compared with the results in Fig. 2 for the one-boson 
truncation. To see the comparison more easily, the one and two-boson contri- 
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Fig. 6. Bosonic structure functions fB±{y), as defined in Eq. (15) of the text, for 
the dressed-fermion state with mass M = fiQ. The wave functions are computed 
from a truncation to two constituent bosons. The longitudinal resolution is varied 
from X = 30 to ii' = 70, with the transverse resolution fixed at = 30. The bare 
coupling is g' = 2. For (a) and (b) the PV masses are mi = fJ,i = IOOO//0, and for 
(c) and (d) they are mi = 10000/io and fii = 100/io- In (a) the amplitude of /b+ 
decreases with increasing K, and in (b) the amplitude of fs- increases. 



butions in each case are plotted separately in Figs. 15 and 16. The differences 
in the results, between the one-boson and two-boson truncations, reflects not 
only the effects of including the two-boson contributions to the kernel in the 
one-boson equation, but also the differences in the bare fermion mass and the 
bare coupling. In the case of unequal PV masses, as in Figs. 15(c) and (d) 
and 16(c) and (d), the perturbative calculation significantly overestimates the 
two-boson contribution. 



In principle, these solutions can be used to estimate the optimal PV mass for 

the one-boson truncation. As the PV masses are increased, we expect the PV 
contributions to the one-boson Fock sector to decrease, but the magnitude 
of the total contribution from the two-boson sector to increase. The latter 
measures the truncation error for a one-boson calculation that does not include 
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Fig. 7. The bare mass niQ of the constituent fermion (fihed circles), the bare couphng 
g (triangles), and the bare- fermion amplitude zq (squares) as functions of the PV 
fermion mass mi . The ratio /Ui / mi of PV boson mass to fermion mass is fixed at 
0.01. The dressed mass is M = /j^q, the radius is i? = O.01//Lto, and the numerical 
resolutions are if = 50 and TV = 30. 

two-boson states even perturbatively. (To estimate the optimal PV mass for a 
two-boson truncation, we would need to estimate the three-boson contribution, 
which would require a much larger calculation.) These two measures, the PV 
contribution to the probability of the one-boson sector, ei, and the ratio of 
the total probability of the two-boson sector to that of the one-boson sector, 
62, are plotted in Fig. 17, along with 62^"^*, a perturbative estimate of 62- We 
compute Ci as the ratio of the PV contribution to the physical contribution, 
each of which is individually divergent but regulated by a transverse cutoff at 
qj_ = (3//.i)^. The fermion mass is held fixed but very large. The optimal PV 
boson mass would be chosen to make the two errors equal; this apparently 
corresponds to a mass on the order of 250//o- 



5 Discussion 



We have solved for the dressed fermion state in Yukawa theory using Pauli- 
Villars regularization and a truncation to no more than one fermion and two 
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Fig. 8. Same as Fig. 7 but for the axial coupling (filled circles), the anomalous 
moment k (squares), the average number of bosons ub (upward triangles), and the 
average momentum fraction (y)^ (downward triangles). 





a 



(b) 



Fig. 9. Bosonic structure functions /b± for the dressed-fermion state with mass 
M = ^0 a-nd radius R = O.Ol/^o- The different curves correspond to increasing 
values of mi. The ratio //i/mi is fixed at 0.01. The numerical resolutions are K = 50 
and N = 30. 

bosons in the Fock expansion. The solution yields the bare coupling g, the 
bare mass mo of the constituent fermion, and the Fock-state wave functions 

as functions of the PV masses fii and mi and the dressed mass M and radius 
R. A limited region of the parameter space has been explored; however, a wide 
range of PV mass values was studied in order to observe the limiting behavior 
as these masses approached infinity. We also studied the effects of truncation 
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(c) (d) 
Fig. 11. Same as Fig. 4 but for the two-boson truncation. 

by comparing the solutions obtained previously with no more than one boson 
in the basis. 



From the wave functions, we calculated various properties of the dressed 
fermion. These are shown in Figs. 12 and 13. The one and two-boson trun- 
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(b) 

Fig. 12. The bare fermion mass mo (circles), the Yukawa coupling g (triangles), and 
the bare fermion amphtude zq (squares) as functions of the radius R. For (a) the PV 
masses are mi = = 2000/io, and for (b) they are mi = SOOOO^o ^ind fii = 500/Uo. 
The filled symbols correspond to the two-boson truncation, and the open symbols to 
the one-boson truncation. The dressed-fermion mass is M = /xq, and the resolutions 
are K = 50 and = 30. For the one-boson truncation, the two-boson contribution 
is computed perturbatively. 
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(b) 

Fig. 13. Same as Fig. 12 but for the axial coupling qa (circles), the anomalous mo- 
ment K (squares), the average number of bosons ub (upward triangles), and the av- 
erage momentum fraction [y) b (downward triangles) . The filled symbols correspond 
to the two-boson truncation, and the open symbols to the one-boson truncation. 
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Fig. 14. Same as Fig. 2 but for the wave functions computed from a truncation to 
two constituent bosons instead of one, with the bare couphng g and bare-fermion 
mass mo adjusted to maintain M = /xq and R = 0.01//io- The resolutions are 
K = 50&ndN = 30. 



cations agree in the weak-coupling, large-radius limit, as expected, but show 
significant differences for stronger coupling, even though a two-boson contri- 
bution is inchidcd pcrturbativcly in the case of the one-boson truncation. The 
presence of the two-boson intermediate states in the kernel of the effective 
integral equation does have important effects. 



The use of light-cone quantization, Pauli-Villars regularization, and carefully 
crafted quadrature schemes can produce accurate solutions for bound states in 
quantum field theory. These techniques should be applicable to more interest- 
ing situations, such as the dressed electron in quantum electrodynamics [12] 
and two-fermion bound states in Yukawa theory and QED. As applications to 
a gauge theory become better understood, we can hope to develop methods 
sufficiently robust to solve for the bound states of quantum chromodynamics. 
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Fig. 15. One-boson contributions to the structure functions fB±{y) as computed 
in the one-boson (dashed) and two-boson (sohd) truncations. The dressed-fermion 
mass \s M = fiQ, and the radius is i? = 0.01///0- For (a) and (b), the PV masses are 
fni = A*i = 2000/xo, and for (c) and (d), they are mi = 50000/io and ni = SOO/xq. 
The two-boson truncation is calculated with resolutions K = hQ and N = 30. 
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Appendix A: Interaction kernels 



The self-energy in the reduced, coupled system Eq. (26) is 
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Fig. 16. Same as Fig. 15, but for the two-boson contribution. In the case of 
one-boson truncation, this is computed perturbatively. 
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Fig. 17. Measures of error in the truncations used as functions of the PV boson 
mass Hi. The quantity ei is a measure of the PV contribution to the probabihty 
of the one-boson sector, and 62 is a measure of the truncation error, defined as 
the ratio of the total probabihty of the two-boson sector to that of the one-boson 
sector. The latter is computed both in the two-boson truncation and perturba- 
tively in the one-boson truncation. The PV fermion mass is fixed at a large value, 
nil = 50000;Uo, and the coupling is fixed at g = 2. The bare fermion mass niQ 
is allowed to vary, to maintain the constraint of = /Xq on the dressed-fermion 
mass. The dressed-fermion radius is allowed to vary. The resolutions are if = 50 
and N = 30. 
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as well as a shift in the azimuthal angle 0' 
can be done analytically; they are 



b' . The angular integrals 



28 



2-K 

I 



277 D + Fcos(f)' ~ v^D^^TF' 



2tt 





27r 



d(p' COS0' - (i ^ 1 

27r + F COS 0' ~ F I ~ ^/W^T^ ) ' 



cos 



2 J,/ 



27r D + Fcos^' 



F2 1 VF*' - F2 y ■ 



(35) 
(36) 

(37) 



The integrals in the self-energy lija also can be done analytically. After a 

— * 

change of variables to ^ = q'~^/q'^ — y'/{l — y) and k± — q'^-\- ^q±, we obtain 
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and /o, -^1, and J defined in Eqs. (23) and (24). For J we have /iQj(-M?) = 
—MjIo(—Mj). For 7o and h, integration over k± yields 
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Integration over z obtains 
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with = {ml -1^1 + M])/{2M]) and Xji^b = y^.^ + /i^/Mf 

In practice one needs to use an alternate form when M| is small, as happens 
when q\ is zero and y is near 1. This alternate form is obtained as an expansion 
in powers of M'j of the integrand in (42), written as 

} InfM^) 

L„i,6 = / dzz"" ln[zml + (1 - z)iJ,l + Mh{l - z)] (45) 

J ■' n-\-l 



The dependence on In [M'^j cancels in the sums in (40) and (41). When 
Mjz{l — z)/ {zm^ + (1 ~ ^)a*6) is of order 0.01 or smaller at z = 1/2, two- 
term expansions of the first term in (45) in powers of M? are sufficient as 
replacements for (43) and (44). 

There is also a special form needed when M| is large, STich as when z is 
near zero. When M| is greater than lOm^ or lO/x^, expansions of the analytic 
expressions (43) and (44) to order /J-l/M^ and \a{Mj/ ^q)ixI/M'^ are used. 



Appendix B: Quadrature schemes 



To solve the integral equations (26), we convert them to a matrix equation via 
quadrature in y' and q'1 and then diagonalize the matrix. For integrals with 
an upper limit of y' = 1 — y, the form of the integrand at 1 — y is obtained by 
explicitly taking the limit; the PV counterterms insure that this limit is finite. 
The quadrature weight is reduced by 1/2 at that point, to take into account 
the edge effect. 

A particularly useful set of quadrature schemes is based on Gauss-Legendre 
quadrature combined with variable transformations to allocate quadrature 
points to important regions and to reduce the q'l integral to a finite range. 
The transformations for the y' integral are 



y'{t) = t\l + dt)/[l + d- {3 + Ad)t + (3 + Qd)f - Mt^ + 2di^], (46) 
t{u)^2u-l. (47) 

The new variable u ranges between -1 and 1, which is the nominal range for 
Gauss-Legendre quadrature. The transformation from y' to t is motivated by 
the need for an accurate approximation to the integral J, defined in (24). 
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This integral is largely determined by contributions near the endpoints when- 
ever the PV masses are large. The transformation y'{t) places many of the 
quadrature points near and 1. It was found empirically, beginning with a 
transformation constructed to compute the integral /[ln(y + eo) —Iniij + ei)\dy 
exactly from the quadrature formula, with eo and ei small. The final form is re- 
stricted to be symmetric under the transformation t ^ {1 — t). The parameter 
d is chosen such that y' ~ O.Olt^ for small t. 



For the transverse integral, the transformation is 

{bya^y-^ - 1 



with V in the range to 1. Only the positive Gauss-Legendre quadrature points 
of an odd order are used for v between -1 and 1, so that v = 0, and therefore 
q'j_ = 0, is always a quadrature point. The points in the negative half of the 
range are discarded. To maintain the accuracy of the underlying sum, the 
weight of the point at zero is reduced by a factor of 1/2. The transformation 
from q'l to v is motivated by its ability to obtain an exact result for the 
integral /[l/(a^ -|- g^) — l/(&^ + q'^)]dq'^. For use in the integral equation (26), 
the parameters a and b are chosen to be the smallest and largest mass scales, 
i.e. a = fiQ and b = uii. 
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